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Abstract 

We summarise the recent theoretical progress in few-body descriptions of the ttNN 
system. Previous descriptions, both three- and four-dimensional, are shown to possess 
serious theoretical inconsistencies. We illustrate how three-dimensional approaches 
suffer from renormalisation problems, and how four-dimensional descriptions contain 
both overcounting and undercounting of diagrams. We then show how such theoretical 
problems have been recently overcome, leading to new practical few-body equations 
for the ttNN system. 

1. Introduction 

In the absence of three-body forces, Faddeev equations provide the theoretically exact 
way to describe three-body systems. It is this fact that has enabled models of the three- 
nucleon system where the effect of the missing three-body forces can be accurately studied. 
It may at first seem that the system consisting of one pion and two nucleons can be similarly 
described in an accurate way with Faddeev equations. This, however, is not the case. The 
problem of course, is that a pion can get absorbed by one nucleon and then emitted by the 
other nucleon, thus making the standard three-body description inappropriate. Indeed any 
number of pions can get created and absorbed by nucleons, and it is clear that field theory 
must be used to describe the irNN system, rather than the standard three-body theory 
of quantum mechanics. The problem of formulating a few-body description for the ttNN 
system, that is the analogue of the Faddeev description for three nucleons, has by now a long 
history. |l| Yet until very recently, both three- and four- dimensional formulations have had 
serious theoretical inconsistencies. Here we would like to summarise these inconsistencies, 
and to describe the recently developed few-body descriptions that appear to overcome all 
these theoretical problems. In particular, we present the new nNN equations in a form that 
is easily compared with previous works, while at the same time being especially convenient 
for numerical solution. 



2. Three-dimensional formulation 



Quantum field theory requires the use of four dimensions to assure manifestly covariant 
descriptions. However, solving four-dimensional equations presents an enormous numerical 
task, and consequently three-dimensional formulations of the irNN system are most desirable 
and indeed have been the most numerous. 

Until recently, perhaps the most sophisticated few-body description of the irNN system 
has been the "unitary NN — itNN model" 0-0. This is a field-theoretic model based on 

1 Invited talk, Sixth International Symposium on Meson-Nucleon Physics and the Structure of the Nucleon 
Blaubeuren/Tuebingen, Germany, 10-14 July 1995. 
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time-ordered perturbation theory, it takes into account pion absorption, and has the desirable 
properties of two- and three-body unitarity. The essential feature of the model is that it 
describes all the processes nd -> nd, nd -> ttNN, NN -> vrAW, NN -> vrd, and AW -> AW 
all within the one set of coupled equations. 

The derivation of the unitary NN—irNN equations is based upon truncating Hilbert space 
to states of at most one pion. In practise, this means retaining all diagrams contributing 
to subsystem nN and AW potentials, but neglecting all other diagrams having two or more 
pions in an intermediate state. Many calculations have been performed with the NN— irNN 
model fl6|-[n|. In general, one can say that the model can account for an extensive amount 



of data, albeit only in a qualitative way. 
2.1 The renormalisation problem 

Despite the modest successes of the unitary NN — ttNN model, it has become clear that 



the model itself has a serious theoretical inconsistency Hl2| , [13fl . The origin of the problem lies 
in the truncation of Hilbert space used to derive the NN — nNN equations. This truncation 
has serious consequences for the renormalisation of both the two-nucleon propagator and 
the ttNN vertex. In Fig. ffl(a) we show the irN nucleon pole diagram where the intermediate 
state nucleon is dressed by one-pion loops; however, the initial and final state nucleons do 
not include dressing since two-pion states are neglected in the truncation. Since close to the 
nucleon pole the dressed one-nucleon propagator is of the form g(E) ~ Z/(E — m), where 
Z is the residue at the pole, Fig. 0(a) illustrates how each irNN vertex f(E) gets effectively 
renormalised by a factor of Z 1 ^ 2 . Thus f W NN — Z x l 2 f{m) is essentially the irNN coupling 
constant, and this fact is used to fix the strength parameter in the form factor f(E). With 
all other parameters of f(E) fixed to reproduce experimental nN phase shifts, this form 
factor then enters the unitary NN— irNN equations as an input. As illustrated in Fig. |l](b), 
when the AW one pion exchange (OPE) amplitude is calculated in the unitary AW — 7rAW 
model, the initial and final nucleons are dressed by pions and consequently each external 
nucleon obtains a renormalisation factor of Z x l 2 . The first renormalisation problem is the 
fact that Z ^ Z. This arises because two nucleons cannot be dressed at the same time in 
the truncated Hilbert space; thus, each nucleon in a two-nucleon state cannot obtain its full 
dressing. This, however, may not be such a serious problem since, in practice, the difference 
between Z and Z turns out to be quite small. The serious problem, instead, is the size of 
the effective 7rAW coupling constant in the AW — irNN equations. Taking Z w Z, Fig. |l](b) 
illustrates that each vertex gets renormalised by a factor of Z, so that the effective nNN 
coupling constant here becomes Zf(m); this is a factor Z 1 ^ 2 times the physical coupling 
constant. With Z being typically between 0.6 and 0.8, we come to the disturbing conclusion 
that the effective nNN coupling constant in the AW — 7rAW equations is smaller than the 
one used in constructing the wN input. This observation helps explain why one typically 
obtains much too small pp -^n + d cross sections using this model [|7|-[|T0|]. 



(a) X Z / ( b ) Z 
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Figure 1: Allowed dressing in the unitary AW— irNN model, with associated Z renormali- 
sation factors, (a) tcN nucleon pole graph, (b) AW OPE graph. 
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2.2 The irNN convolution equations 



Here we describe a new formulation of the ttNN problem where unitary equations are 
obtained without having to truncate the Hilbert space to some maximum number of pions. 
Consequently, all possible dressings of one-particle propagators and vertices are retained in 
our model. The essential technique that enables the calculation of all such dressings is a 
novel use of convolution integrals. In this way we overcome the renormalisation problems 
discussed above. 



As an explicit derivation of the new irNN equations can be found in Ref . ||14|| , here we 
prefer to simply state the final equations and to describe their essential features. The new 
ttNN equations can be expressed in many different forms, all of which are equivalent. The 
form we shall choose here is the one that most closely resembles the unitary NN — irNN 
equations as given by Afnan and Blankleider|5|] (AB). Choosing this AB form has two es- 
sential advantages: firstly, we are able to directly compare the differences between our ttNN 
convolution equations and the unitary AW — irNN equations, and secondly, this form is 
ideal for numerical solution, especially since advantage may be taken of existing codes for 
the unitary NN — irNN equations where essentially the AB form has been used. We note, 
however, that for the sake of easy comparison, we give the equations here only for the case 
of distinguishable nucleons. Including proper antisymmetry of the nucleons is a technical 
formality which, for the AB form, will be presented elsewhere (in this regard, note that the 





derivation of the four-dimensional irNN equations presented in Ref. pj], assumes identical 
nucleons from the beginning). 

The irNN convolution equations may be expressed as a set of coupled equations for the 
reactions NN — > NN, NN ird, ird — > NN, and ird — > ird using the following (4 x 4) matrix 
form: 

f rr> rh \ f t/ t? \ ( / ri n \ / m rh \ ^ 

(1) 

Before explaining the symbols in this equation, let us define what we mean by a product 
of two symbols. For any two quantities B and A, describing processes m k and k — ► n, 
respectively, we define the product symbol AB to mean the the integral 

AB = I dp'l . . . dpi A(p[ ...p' n , p'[ . . . pi; E) B(p'i . . . p' fc \ Pl . . . p m ; E) (2) 

where pi is the three-momentum of particle i and E is the total energy. Although, momentum 
conserving ^-functions are assumed to be contained in both A and B, it is easy to see that 
such 5-functions can be factored out without affecting the symbolic equations. In Eq. ([[]) 
the unknown quantities are T NN , together with T aN , T Nf3 , and T a/3 which are elements of the 
matrices T_ N , T_ N , and T, respectively (here indices a and (3 take on values 1, 2, and 3). The 
physical amplitudes for NN—>NN, NN^ird, ird^NN, and nd—>-7Td, are then given by 

Xnn — Tnn ; X dN = \l/rfT 3A r ; X Nd = T N3 ^ d ; X dd = ^/ d T 33 ^ d , (3) 

respectively, where ty d is the deuteron wave function in the presence of a spectator pion. 
On the r.h.s. of Eq. (|l|) G is the fully dressed irNN propagator, G^n is the fully dressed 
NN propagator, I is a 3 x 3 matrix whose elements are 8 a p = 1 — 5 a p, and Vnn is the 
dressed one-pion exchange potential given by Vnn — Y%j=i F&jGoFj where F« and Fi are 
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Figure 2: The vertices Fi and Fx- The dark circles represent all possible intermediate states 
consistent with the requirement that F\ and Fx be amplitudes with chopped external legs. 
Vertices F 2 and F 2 are obtained by interchanging 1 and 2. 

fully dressed irNN vertices in the two-nucleon sector as illustrated in Fig. |2|. Finally we have 

the matrices , n n _ N 

( wi w% \ 

Z=(E?=iM>) ; £=(e?=i^/j) ; «>§ ^2° o (4) 

V o o W3 ° ; 

where the w° (a = 1 . . . 5) are the disconnected iVTV-irreducible amplitudes for irNN —*7rNN, 
to be discussed shortly. 

By form, Eq. ([[[) is very similar to the unitary NN — irNN equations as given in Eq. 
(59) of AB. However, the essential feature of Eq. (|l|) that distinguishes it from the unitary 
NN— irNN equations, is that all input quantities in Eq. ([T]) are fully dressed. In this way the 
renormalisation problems of the NN — irNN equations have been overcome. However this 
would only be a formal solution to the renormalisation problem if it were not for the fact that 
all the necessary dressings can be calculated exactly using convolution integrals. That this 
is so follows from Ref. jn| where we showed that any disconnected Green function is equal 
to the convolution of all its disconnected parts; thus for example, the dressed two-nucleon 
propagator Gnn is expressed in terms of the dressed one-nucleon propagators gx and g 2 as 

1 r°° 

G NN (E) = - — / dzgx(E - z)g 2 {z) (5) 

where, for the sake of simplicity, we have set the momenta of the nucleons to zero. To further 
save on notation, we introduce the shorthand Gnn = 9i®92 to mean the convolution integral 
of Eq. (|). Giving labels 1 and 2 to the two nucleons, and label 3 to the pion, in the same 
way we have that the fully dressed irNN propagator Go is given by the double convolution 

G = g\ ® g 2 ® g s . 

To see how the amplitudes u>° are calculated, we first define the amplitudes w a to be 
the disconnected irNN — > irNN amplitudes, illustrated in Fig. each corresponding to a 
different type of disconnectedness, and containing all possible contributing diagrams. It is 
just because the w a contain all possible contributions that one can express them through 
convolution integrals as 

Wx = h g> g 2 ; w 2 = t 2 ®gx ; w 3 = t 3 ® # 3 ; w 4 = f x ® f 2 ; w 5 = / 2 ® f x (6) 



1 9 ~— _ 1 
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Figure 3: The amplitudes w a . The dark circles represent all possible intermediate states 
consistent with the requirement that the w a be amplitudes with chopped external legs. Am- 
plitudes w 2 and w 5 are obtained by interchanging the two nucleons in wx and W4 respectively. 
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where the "tilde" denotes a Green function quantity consisting of the corresponding ampli- 
tude with additional initial and final-state propagators; thus, for example, w a = Gow a Go, 
and t\ = 5 f 7 rAr 1 ti5' 7 rAr 1 where t\ is the t-matrix and g^Nj the dressed propagator for scattering 
of a pion off nucleon 1. As we have shown in Ref.flnj, the convolution integrals effectively 
sum over all the relative time orderings of one subamplitude of a disconnected diagram with 
respect to another. Similarly, the vertices F\ and F 2 are also expressed via convolutions as 

Fi = h®g* ; F 2 = f 2 ® 9l . (7) 

Once the w a are calculated, we may then write them as w a = u>° + where is the part 
of w a that is two-nucleon reducible, while W® is two-nucleon irreducible. Since we consider 
all possible contributions, it is clear that (i = 1, 2) wf = FjG^Fi, = 0, wj[ = FiG^F 2 , 
and w§ = F 2 GnnFi. In this way, all the essential input to Eq. (|l|) has been specified. 

We may finally note a second major difference between Eq. ([I]) and the unitary NN—ttNN 
equations. The input matrix w° in Eq. (|l|) has off-diagonal elements, while the corresponding 
matrix for the NN — irNN equations is diagonal. Recalling that the amplitudes of w° are 
two-nucleon irreducible, we can see from Fig. that the off-diagonal elements u^j and 



correspond to what has been called the Jennings terms. As pointed out by Jennings Hl6|, these 
terms may be important for the understanding of ird scattering. In our case, the Jennings 
terms are also fully dressed, and form an essential part of the convolution equations. Indeed, 
since our NN propagator Gnn is fully dressed, it also contains two-pion states coming from 
intermediate Jennings-like terms. It is then necessary to retain and in the convolution 
equations because they combine with Gnn in just the right way to guarantee three-body 
unitarity. 

We recall, that the only approximation made in deriving the convolution equations of Eq. 
([I]) is the neglect of all connected 7riV7V-irreducible diagrams for the ttNN —>7tNN process ||T4|. 
Yet it is very easy to include some types of connected contributions. One such contribution 
would involve intermediate state potentials that are 7riV7V-irreducible. Then Eq. ([!]) 
would be modified simply by replacing Vnn with Vnn + Vnn- This observation suggests that a 
way to include heavy meson exchange into our NN potential would be as a phenomenological 
model for V^f. 

3. Four-dimensional formulation 

Although three-dimensional equations may be easier to solve than those in four- dimensions, 
there are important reasons why the formulation of four- dimensional equations is necessary. 
Firstly, such equations are based on relativistic quantum field theory, and retain the fun- 
damental property of off-shell covariance. Secondly, having the correct four- dimensional 
equations, one can then do a three-dimensional reduction using one of the well-known re- 
duction schemes. We may also add, that with the ever increasing power of computers, the 
numerical solution of four-dimensional equations becomes ever more feasible. 

The first attempts to formulate few-body equations using relativistic quantum field theory 
were made already in the early 1960's |T7|]-[[L9(|. Both such general formulations and ones 



more specific to the nNN system have been pursued until the present time ||20||- ||23|| . Yet 
as in the three-dimensional case, all these attempts have had theoretical inconsistencies. 
In particular, all previous attempts have contained either overcounting or undercounting of 
Feynman diagrams. 
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Figure 4: Example of overcounting in NN —> nd. (a) The NN —> nd Feynman diagram 
where dark circles represent all possible contributions, (b) One of the contributions included 
in (a), (c) Another way of drawing diagram (b) showing how overcounting arises. 

3.1 Overcounting and undercounting problems 

Perhaps the easiest way to illustrate the overcounting problem in the nNN system is 
with an example. Consider the "triangle" diagram of Fig. |](a) for the process NN —> nd, 
where the dark circles represent the full nN ^nN amplitude, the dressed nNN vertex, and 
the dressed deuteron vertex. If one were to calculate this diagram in four dimensions, as is, 
using covariant forms for the off-shell nN t-matrix, nNN vertex, and the deuteron vertex, 
then one would have the mistake of overcounting of diagrams. This is illustrated in Fig. |](b) 
where we consider just the crossed-pion graph contribution to the input nN t-matrix. As 
these are Feynman graphs, there is no meaning associated with the slope of the lines, and 
one could just as well have drawn Fig. f|(c). However Fig. f|(c) clearly illustrates that this 
contribution corresponds to the dressing of the already fully dressed deuteron vertex. 

This type of overcounting arises in four- dimensional approaches whenever one tries to 
formulate multiple-scattering graphs in terms of fully dressed vertices and full amplitudes 
for all subprocesses. In once-off cases, like that of Fig. [|(a) , one can easily fix the overcounting 
problem by making a necessary subtraction (here one would subtract the graph of Fig. §(b) 
from the calculation of Fig. |](a)). However, the way to solve the overcounting problem for the 
case of coupled integral equations is highly non-trivial as an infinite number of overcounted 
contributions are involved. 

In a similar way, let us illustrate how undercounting arises in the covariant nNN problem. 
As in the three-dimensional formulation, one neglects three-body forces also in the four- 
dimensional case. Only in this way can one obtain few-body equations where (in the cm.) 
no more than two independent momenta are involved. However, one does need to be very 
careful about neglecting three-body forces in the four-dimensional theory. Consider, for 
example the Feynman diagram of Fig. ||(a). This is a graph for the process nNN —> nNN 
that is both connected and n NN- irreducible. It therefore corresponds precisely to what is 
meant by a three-body force. However, neglecting this contribution from a few-body theory 
of the nNN system would be a bad mistake. This is illustrated in Fig. |5](b) where we allow 
the graph of Fig. ^(a) to couple to the NN channel. Again no meaning can be attached 
to the slope of the propagator lines, and we can equally well draw this diagram as in Fig. 
D(c). This, however, reveals that the three-body force of the nNN — > nNN process has now 

(a) T (b) 



Figure 5: Example of undercounting in NN — > nNN. (a) A nNN — > nNN graph that has 
usually been neglected since it corresponds to a three-body force, (b) The coupling of the 
graph in (a) to the NN channel, (c) Another way of drawing diagram (b) reveals a two-body 
process. 
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become a two-body rescattering contribution in the NN — > irNN process. Thus neglecting 
the three-body force of Fig. ||(a) would lead to an undercounting of important two-body 
contributions. 



3.2 Four- dimensional ttNN equations 

In a recent paper, we have solved both the overcounting and undercounting problems 
in the formulation of few-body equations in field theory |24] . The few-body equations for 



the ttNN system then follow as a particular case. The method used to derive the equations 
involves the classification of Feynman diagrams according to their irreducibility. The over- 
counting problem is handled by a procedure where, in formally identical cases like that of 
Figs. |](b) and (c), one of the two right-most vertices is "pulled out" further to the right. 
The undercounting of diagrams is handled simply by retaining all three-body forces until the 
end of the derivation where the ones that did not lead to two-body interactions are safely 
neglected. It is gratifying that Phillips and Afnan|25| have recently confirmed our equations 



using a modified version of Taylor's original classification of diagram scheme (T^, 26fl. 



With three-body forces neglected as described, one might still find it useful to retain, as 
in the three-dimensional case, the NN—>NN 7riV7V-irreducible potential vj^, as well as the 
simultaneously NN- and 7riV7V-irreducible connected NN — > irNN amplitude F n . However, 
let us at first consider the simplest case where these contributions are neglected (they are in 
fact completely absent in the usual case of a (fnpip interaction). 

In this case the four-dimensional nNN equations can be written as Eq. (|l|), but with 
the following modifications: (1) The product of two quantities A and B is defined as in Eq. 
(0), but now with all momenta and integrations being four- dimensional (in particular the 
replacement dp'- — > d 4 pi/(2ir) 4 needs to be made), (2) all convolutions of Green functions 
are replaced by usual products, (3) the matrix w° is now diagonal (i.e. with = w® = 0), 
and (4) the following replacements are made, 

E^F-BL ; F^F-& ; V NN - V NN - A (8) 

where the terms Bl_ and A are subtraction terms that exactly compensate all the overcounting 
due to the use of full off-shell amplitudes and fully dressed vertices in the coupled scattering 
equations. They are defined with the help of Fig. || as follows: 

A = W n7T + Wl N + W NN + X + Y' - B'G B' (9) 
(a) A (b) X (c) 



(d) >< (e) (f) 

i 1 2 ' 1 2 

Figure 6: The subtraction terms in the four-dimensional ttNN equation: (a) W n7T , (b) 
WnN, (c) Wnn, (d) X, (e) Y, and (f) B. The dark circles represent the following two-body 
amplitudes: (a) full tttt t-matrix, (b) one-nucleon irreducible ttN t-matrix, and (c) full NN 
t-matrix minus the NN one-pion-exchange potential. 



x , — - — 2 
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where the dashed quantities are the sums W^ N = W n N + PW^nP, Y' = Y + PYP, and 
B' = B + PBP, P being the nucleon label exchange operator. &_ is a column matrix with 
each element given by B' . 

In the more general case where vj^ and F n are retained, it turns out that only the 
subtraction terms need be modified. In particular, we need to do the replacements A — > 
A - VS - FAFi + F 2 ) + F 2 ^F X F 2 - (F 1 + F 2 )F 7r - F 1 F 2 F 27T and B' -> B' - F v , where F 2w 
is the connected NN- and ttNN- irreducible amplitude for NN — > ttttNN. 



> Permanent address: Mathematical Institute of Georgian Academy of Sciences, 
Z. Rukhadze 1, 380093 Tbilisi, Georgia. 
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